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Использование циклических сдвигов 

для ускоренного вычисления циклических 
сверток длины равной произведению 
взаимно простых чисел 


п-1 

В работе рассматривается метод вычисления циклических сверток длины М = | [7 › где (т,,т,) =]: 
1=0 

ТЕ 7 . Рассматривается новое блочно-циклическое представление свертки за счет использования 


циклических сдвигов в каждом блоке. При таком подходе китайская теорема об остатках (КТО) не 
используется. Метод является эффективным за счет использования несложных пред- и поствычислений 
с применением циклических сдвигов. Он позволяет распараллелить вычисления. Получены априорные 
оценки сложности по количеству операций умножения. Приведены примеры вычислений сверток. 


Введение 


Циклические свертки широко используются при построении цифровых фильтров, 
в задачах ассиметричной криптографии и т.д. Как известно, выбор метода зависит от 
области его эффективного использования. Циклическая свертка используется при 
реализации быстрой операции умножения с использованием быстрых преобразований 
Фурье, Уолша, Хаара [1]. Операция умножения занимает большую часть вычислитель- 
ного времени при реализации операций ассиметричной криптографии, в частности в 
задачах распространения секретного ключа, шифрования информации, выработки и 
верификации электронно-цифровой подписи и т.д. 

В работе [2] Д. Питасси предложил метод вычисления циклической свертки 
длины М№М=2”, усовершенствованный впоследствии Девисом [3]. Данный метод 
приводит операцию вычисления циклической свертки длины 2" к вычислению трех 
циклических сверток половинной длины 2"". На каждой итерации количество 
сверток половинной длины увеличивается в три раза. Итерации продолжаются до 
получения сверток длины 2. В работе [4], следуя Питасси, показано, что вычисление 
циклической свертки длины №М=2.К, К — нечетное, сводится к вычислению двух свер- 
ток длины КА, что позволяет реализовать операцию умножения двух К -разрядных 
чисел, К — нечетное. Метод, рассматриваемый в данной работе, является продолже- 
нием работы [4]. В работе [5] приведен метод умножения больших чисел, основанный 
на циклической свертке длиной, равной степени двойки с использованием прео- 
бразования Уолша. В работе [6] приведены оптимальные алгоритмы по количеству 
операций умножения, сложения и вычитания для вычисления сверток длины 2, 3, 4, 
5, 7, 8, 9. Р.К. Агарвал и Ч.С. Баррас предложили отображать одномерный массив в 
многомерный, используя удвоение числа точек [7]. В работе Винограда [8] для 
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вычисления длинных сверток применяются билинейные формы. В работе [9] Р.К. Агар- 
вал и Дж.У. Кули применяют китайскую теорему об остатках (КТО) для преобра- 
зования одномерной циклической свертки в многомерную, которая циклична по всем 
направлениям. Упомянутые выше работы направлены на уменьшение числа однословных 
умножений. Сегодня архитектура современного процессора такова, что операции одно- 
словного умножения и сложения выполняются за одинаковое число тактов. Поэтому 
дальнейшие работы должны учитывать этот факт и также быть направлены на 
оптимизацию количества таких однословных операций, как сложение, вычитание, сдвиги. 
Одним из оригинальных подходов для решения таких задач можно назвать комби- 
нацию программного и аппаратного подхода, при котором перед началом каждого 
вычисления в зависимости от длины входных данных прошиваются эффективные для 
данных длин подпрограммы [10]. 

Для вычисления свертки длины Л стандартным методом необходимо выполнить 


2 > > 
№ операций умножения, что соответствует умножению квадратной матрицы на 
вектор длины №. Продолжая идею Питасси [2], [3], операцию вычисления циклической 


М№МЬ-—1 
свертки длины № =15 (и, = а И = 0, М-—1) можно свести к вычислению 
т=0 


трёх сверток длины № =5 (рис. 1). 
Хх | № № № №1 № Ха М № 1№0 № № М ФХз 


У № № № № №№ № № 2 м № М м М 


| 

Рисунок 1 — Вычисления свертки длины № =К.-и=3.5, где А, г — нечетные 
№ ... Х. №; ... № Хо ... Хз 

Свертки : ` :,: "1, : *  : повторяются по три раза. 
Х1> ... № х› ... №; Х7 ... №0 


Данный подход позволяет значительно уменьшить число операций умножения, 
сложения, вычитания и т.д. 
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Далее будет показано, что вычисление любой циклической свертки длины 
№М=К.и, (К,г)=1, К<г, можно представить в виде К блоков, повторяющихся К 


раз. Такой метод позволяет свести задачу вычисления свертки длины М=К.г к 
задаче вычисления свертки длины КА. Далее будем использовать выражение 
«секция» вместо «блок», что созвучно с методом парисекции, рассмотренным 
Девисом [3]. 

Новизной предложенного метода является использование циклических сдвигов 
в каждом блоке для представления свертки в блочно-циклическом виде без ис- 
пользования китайской теоремы об остатках. Преимуществом данного метода является 
то, что не используются предвычисления, на каждой итерации достаточно знать 


Р (М=К.г, Е<к, а=Р-Ё, (р.4), =1), не используются дискретные преобразо- 


вания, метод разбиения подобный для любых длин сверток. 


Операторы и обозначения 
Перед тем как рассмотреть предложенный метод, введем следующие операторы 
и обозначения. Последовательности Ху и У, длины М№ представляются векторами- 


столбцами: 


Хм-1 Ум 
Циклическую свертку двух сигналов Ху и У, обозначим оператором ®: 
и м1 


К =Х, ®У,, К =| | пе хьУиьл, = 0,М-1, 


т=0 


№ УФ У У У У | № № ХХ М 
У У У У№ м № № ХФ 
Х› | У У У У ИЛИ ух м № Х.. 


% У № У У У № ХХ М № 


Введем следующие операторы: 5, Р, И, О (5-“беесР” (выбрать), Р-“Рат!” 
(часть), О-“Ор” (вверх), О-“Ром’” (вниз)), и обозначим их следующим образом: 


5 Хы = хан» ЛЕ" -1, 1=0, 1, где К = М/к; 
циклический сдвиг элементов вверх — Г, =ИХу, у, =х (в)? = 0, М-1 : 
циклический сдвиг элементов вниз — И, =ДОХу, у, = Хм, К=0,№М-1; 
циклический сдвиг элементов вверх р раз- Г, =И? (х м ), УЕ Жинр), › К = 0, М-1 : 
циклический сдвиг элементов вниз р раз- Г, =? (х ы ), УЕ Хунм-р), К=0,М-1. 
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Проиллюстрируем операторы для длин М=5,12,15: 
Е ЕЯ Гх, | Гх, | 
№ х> Хх. Хх 
РХ, =|х | ИХ; =| х, |› БАХ =| ж | (ХЗ =|х, |› 
р а р № 
[Хз | [№ | [№2 - Е - 
№ 


Х 
Е Е Г] Жо [| 
№ Хх. х; х № 
5-Х =| № |? 55 5 =| % |? ХЕ =| № |? РЖ =, |› в: =| %) |? 
№ Хо Хи х: № 
[№12 [Хз [№14 [4 | [№ _ 
о 


3 


4 


Постановка задачи. Если целые положительные векторы Ху, У, Ку вида 
М | М м 
ры © РЯ @ —- @ @ 57 
Ее) : = У. У2 ‚ ВА =) 72°, где 0<х,у,к<2, М=7’, тогда 
1=0 1=0 


1=0 


выражение К» вида 


1=0 


М№М-—1 | М№М-—1 . М№М-—1 
ВЫХ, =: [2 ни, 
1=0 


где и. = У представляет собой циклическую свертку Ху и У,. 
р+д=т(тоа М) 


Необходимо найти эффективный алгоритм вычисления циклической свертки. 


Свойства последовательностей 


Лемма 1. Пусть последовательность „Лу составляют элементы /=0, М-1, где 


№М=К.и, К<х, (К,г) =1. Если Лу разбить на К секций длины / с использованием опе- 


ратора 5 , то каждая секция 5'.Л у ‚ {= 0, — 1 , содержит индекс, которой делится на г. 

Доказательство. Так как каждая секция является результатом применения 
оператора 5, то соседние элементы внутри каждой секции отличаются на К. 
Нулевая секция 1=0 будет состоять из следующих элементов: 


О, А, 2%,....("-1)-К. 
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Элементы секции 1=1 будут следующими: 
1,1+^,1+22.,...,1+("-.К. 
Их можно выразить следующим образом: 
2=1+р.К, р=0,"-1, (1) 


где г — длина секции. 
Найдем значение р, при котором 2 делится на г без остатка. 


Сделаем замену А=и-@ в формуле: 
2=(1+р.^),=(1+р-("-а)), =(1+р-"-р-4),=(р), -г+(1-р.4),, р=0,М-1. 
Отсюда ВИДНО, ЧТО 2 будет делиться на г только в том случае, когда 


(1-р.4),,, делится на г: (1-р.а) =0, 1=(р-а) =(9"-а), т.е. когда р=а'. 


| 
Число р можно найти, используя расширенный алгоритм Эвклида: 
4: р-г."=а.а'-г.и=1,г>@а. 
Для случая 1> 1 соотношения (1) можно записать в следующем виде: 
2. =(1+(-р). К), =(А+р.®)), =(-2),. (2) 
С учетом того, что число секций (длины и) равно К и порядковый номер в секции 
не может превышать длины секции (1. р) < г, соотношение (2) примет следующий вид: 


д=1+р-К; 2. =(@+р,.К),=@-2),, р. =(-р),1=0,К-1, (3) 
где А — количество секций длины г, р — порядковый номер элемента в секции 
Г со значением 2,, которое делится на г без остатка. 
Окончательно получаем: 


р=4" =(@-Ю”); р=(.р),,1=0,Е-1. (4) 
д=1+р.к; 2,=(-2),,1=0,-1. (5) 


Лемма доказана. 


Соотношение (5) можно представить в виде 2,=Г. ( . 2/"), ‚ 1=0,К-1, 
упрощающем вычисление 2,. 

Лемма | устанавливает зависимость между числом циклических сдвигов р 
секции, номером секции 1, числом секций К и длиной секции г. 

Для примера рассмотрим последовательность „Л; = (0,1,..., 14), М=15, А=З, 
г=5. Вычислим р=а“" ={@-Ю7) =(2—). =3, 2=1+3.3=10. Найдем р,, 2, 


1=0,^—1=0, 2, используя соотношение (3). 


Таблица 1 — Вычисление 5".Л‚, р, =(7-р),, 2. а 1=0,К-1, 


для М=ЕК-и=3.5=15 , где К, г — нечетные 


=. 


5,„Ты р. =(1-р), ВА, 
5%. ‹ = (0,3, 6,9,12) | р. =(0-3). =0 |2. =(0.1 
ОЕ Ч, рее 9 
бл, 5 ре || =210) 55 


| 
> = о 


=.| =. 


—. 
| 
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Рассмотрим .Л.; = (0,1,...,34), М=35, А=5, г=7. Обратный элемент р будет 
равен р=а'=((#-Ю7) (27) =4, ==1+4.5=21. Найдем р,, 2, 1=0,К-1=0,4, 


г 


используя соотношение (3). 


Таблица 2 — Вычисление 5'.Л‚, р, =(2-р), 2 = (1.2), м 
для М=К-и=5.7=35, где К, г — нечетные 
5.7 рр), д=(Г. 2), 
1=0 | 57Л,; = (0, 5,10,15,20,25,30) | рр=(0-4) =0 |25=(0-21), =0 
1=1 | 51Л,; =(Ь 6,1116, 21, 26,31) р =( 


#1=2 | 52/,, = (2,7,12,17,22,27,32) | р. =(2 
=3 | 531, = (3, 8,13, 18,23, 28,33) | р. =(3- 
[=4 | 57; = (4,9,14,19,24,29,34) | р.=(4 


Анализируя табл.1, 2 и лемму 1, приходим к следующему утверждению. 


Лемма 2. Пусть последовательность „Лу составляют элементы /=0, М-1, где 
№М=К.и, К<к, (Кг) =1. Если У» разбить на К секций длиной / с использованием 
оператора 5, то первый элемент секций И”'5 ей у, = 0.6 —1 ‚ содержит индекс 2,, 
которое делится на г без остатка, где р, и 2 определяются следующими 
выражениями: 

р=4"=(("-Ю") ; №=0; дер; р =(ра+р),,1=2,Ж-1. 
Е=1+р-К; 5 =0;4=2;2=(2..+2),1=2,К-1. (6) 


Доказательство. Применяя рекуррентные соотношения к соотношениям (4) и 
(5), приходим к искомым формулам. С учетом леммы 2 табл. | может быть 
представлена в следующем виде. 


Таблица 3 — Вычисление УЛ. ь ИР ‚ р=25 =0, р, = (р +р) 


г? 


2. и 1=ЬА-1 для М=К.-И=3-5=15 , где К, г — нечетные 


5,7 р. =(7- р), ОРЛ, 2. =(1-2), 
0 15°. = (0, 3, 6, 9,12) р = (0-3). =0 (°$5.Л; = (0, 3, 6,9,12) |2 =(0-10). =0 


11| 55 = (1,4, 7,10,13) | р, =(1.3). =3 |055; = (40,13,1,4,7) |2, = (1.10) =10 
2 | $5.Л; = (2,5,8,11,14) |р,=(2.3), =1 |0'55Л; = (5,8,11,14,2) |2, =(2.10),=5 


Из табл. 3 видно, что секции можно строить сразу, опуская циклические сдвиги 
вверх. Для этого достаточно знать индекс первого элемента каждой секции. Индексы 
остальных элементов в каждой секции отличаются на шаг А, равный количеству 
секций при разбиении. 
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Рассмотрим построение секций, начиная с первых элементов секций, на 
примере последовательности „Л. = (0,1,...,27). М№М=28, К=4, г=7. Вычислим 


р=4"=((-Ю”) =(3") =5, 4-р-г-и=3-5-7.2=1, 2=1+р-к=1+5-4=21. 


Найдем 2,1=0,К-1=0, 3, используя формулу (6). 


Таблица 4 — Вычисление лы. ‚ р.= (; . РЕ 2, = (: . бе ‚1=0,А-1, 


для М=К-г=4.7=28, где К — четное, г — нечетное 


7 _ = (1-2) 0"5 Л» Р; =(#-Р), 
вы: О 075.21 25,55,0.13:17). | р=5). = 


1=2 |2,=(2.21),=14 | (*51Л» = (4,18, 22,26,2,6,10) | р, =(2.5), =3 
(3.21, =7 (187.7 = (7,11,15,19,22,27,3) | р, =(3-5), =1 


Таблица 5 — Вычисление секций для М=К.и=9.11=99, 
где К, г - нечетные, р=6, 2=55 


: Е = (1-2), | ы. = (1-2), у 2: =(-2) 
0 5 =0 3 2, =66 6 ге 98 
1 #=55 4 = 22 7 2. =88 
2 2, =11 5 ЕЙ 8 2, =44 


Лемма 3. Пусть последовательность „Л = (0,1,..., М-1), М=К.г, (Кг) =1, 
к <г упорядочена: Р'(Т,) =И”($'(,)), р, = (1. р).,1=0,К-1, где р=(("-Ю”) 


Тогда к первоначальному порядку полученную последовательность Г, переводят 


т 


следующие соотношения: 
р=((-Ю") ; РТ.) =Б®(РАТ,)), р=(-р),,1=0.Ж-Т. 


ВЕРЕ РЕЛЬ (7) 


Доказательство. Оператор И” (ир-вверх) в выражении И” 5"./, циклически 
сдвигает вверх каждую секцию 7 на р, раз. Соответственно, чтобы получить зна- 
чения секции в первоначальном (исходном) состоянии, нужно применить обратный 
оператор ДО (4оим-вниз) для каждой секции р, раз, т.е. циклически сдвинуть вниз 
каждую секцию 1 на такое же р, раз: Р/Т, = О*'(Р/Т,), 1 =0,к -1. После этого перво- 
начальными элементами каждой секции 7 будут элементы со значениями 1, которые 


отстоят друг от друга на г элементов (длина каждой секции равна г). Получили следу- 
ющую последовательность, состоящую из секций, выделенных вертикальными линиями: 


СОТ МЕ 
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Выборка первых элементов, отстоящих друг от друга на ’ шагов, из К секций 
определяется выражением 5.Т) ‚ что соответствует последовательности (0,1,....К-—1). 
Выборка элементов (К, А +1,....2К-1), следующими за первыми в каждой секции, 
определяется выражением 5,Т,. Продолжая выборку по формуле (7), получим 
первоначальную последовательность. 

Лемма доказана. 

Более детально рассмотрим последовательность Л; длиной 15 элементов 


К=3, г=5. Согласно табл. 3 упорядочивание Р,Л.; =И”'5.Л:, р, = (1.3) ,1=0,2, 


г 


дает следующую последовательность, разбитую на 3 секции: 
Т; =(0,3, 6,9,12 | 10,13, 1,4,7 | 5,8,11, 14,2). 


Выражения РТ; =02”(РЛ.), р,= ( : 3) 1=0,2, приводят предыдущую 
последовательность к следующему виду: 


Т; = (0,3, 6,9,12 | 1,4, 7,10,13 | 2, 5, 8, 11,14). 


5 > 


Соотношения Р/.Л,; =53Г., 1=0,4, приводят к первоначальной последова- 
тельности: 

Л; = (0,1,2 |3, 4,5 | 6, 7,8 | 9, 10,11 | 12,13, 14). 

Лемма 4. Циклическую свертку А» длины М=К-и, К=г-1, (К,Г) =1 сигна- 
лов Ху и У, выражения РХ,=И'5'Х,, РУ, =0'5'У,, РВ, =0'5' Ву, 1=0,К-1 
переводят в блочно-циклический вид. 

Доказательство. Циклическая свертка сигналов Ху и У, определяется 

М№М-—1 


выражением Г, — Ух 


т=0 


‘У, Л=О0, М-1. 


т) 


Воспользуемся оператором 5 и найдем секции длиной г для сигнала Ху» 


(преобразования для сигнала У’ не приводятся, т.к. они аналогичны): 


№ м Хе 
У Ху Хо 
ео А 9 || (8) 
Х м2 Хм Хм 
Е мк Е м-н [Хм | 


Элементы каждой секции 1, полученные с использованием оператора 5, 


сдвинем циклически вверх на 1 раз 1 =0,^ —1: 


№ Хх Хх Хх, 
Хь нь х, Х к 
00 тс! 1 1 
И 5,.Ху=| ... 0 5.Ху=9Ч =( зы т ых + 
Х м2 Х+м-2% А+ м-2% Хм 
Хм | Ен Ен х 
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о ы х», 
>44 Х+^ 2+ 
№ О |= 
выю, ИИ 2+2 у 2" > 27+2к 
[2+ № [мк [1 №2 | 
х,, (к-р 
Хы р-н 
Рей _ с ф — —_ 
05,Х у = но, Ар 2 = 
Х 24 Хе 
Е Жк [ Ж-ьни-к 


э 
э 


Хм 


Хм + 


Хм-„-2к 


| Хм 4 


8Т 


(9) 


Из (9) видно, что индексы при х первых (самых верхних) элементах каждой 
секции делятся на г: 


[95 Хи =, =0,К-1. 


г) 


Так как каждая секция построена таким образом, что индекс при х следу- 
ющего элемента в секции отличается от предыдущего на А, то все элементы секций 
можно записать следующим образом: 


[015 Х,],= 


<2;+]-К>м 


2 = 0,"-1, 2, =1-Г, 1=0,К-1. 


(10) 


Вычисление свертки с учетом (10) можно представить в следующем виде 


АЕ 7 


ПаинтК) = 83 


1 
- еек уж-т к), У (ЛК) 
50 /= 


‚т=0, 7-1, 4=0,К-1. 


(0 


Из данного соотношения видно, что в каждой колонке (4-г + т. К) матрицы 


(рис. 2) индексы при х и у каждой строчки ( "7+: К) х отличаются на одно и то же 


число (- д:г-т-К ) „у › Что говорит о свойстве цикличности данного представления. 


Рассмотрим вычисление свертки. 


| | | 
У № Хм Хх | Хм, Хм Хм Е м + 
| г 
У+ и № Хх Хм Хм р В к х, Хи 
| ИХ 
| | | 
| г 
| г 
Ум Хмк Хм № р Хм Хм Хм, м Ам х, 
Ее ОЕ 
У, х, хм Хы № Хм Хх КО о Ян + 
[ [ | 
Ук Хы х, 4 Хнок Хх №0 Хх т № Х>+2к 
| Вы 
| .. г 
| Е 
[ | | 
НИ М Хм х, р Хм-к _ Хм № 2 22 Хх», 
О р р о = О О ат ни 
ри а ВВ ее Ме ОИ пар ее 
[ | | 
Ук-, Хм, Хм Хм-1+к Хм. Хм Хм -21+к г № Хм Хх 
[ | | 
Ум-и+к Хм Хм Хинок | 1| Хм Мм» Хм | м №0 Хх 
| Г... | 
| | | 
| РЕ 
| И 
Ун-г-к Ху--к Хм Хм, [Хм ХМ, Хм [Хм Хм № 
| | 
Г" о 7, Пу | 7, 7. п | [. м Гу | 
Рисунок 2 — Вычисление свертки длины М=А-и, Г=К-1 
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Для доказательства блочности представления достаточно исследовать верхние 
левые элементы секций (блоков). С учетом соотношения (11) (при т=0, 1=0) 


получаем следующие выражения: 


У» =, -1, 4=0,#-1, (12) 


х, "4" 


где й ” У: ) элемент строчки (; Е г), И столбцы (- Ч: го матрицы. 


мат) ТИ 


Откуда видно, что угловые элементы х 1=0,^-1, 9=0, -1, отстоят 


(1-4) у 2 
друг от друга на г строчек и столбцов. Выражение (0 -4)-"), =г. ( —4) в 20, К-1, 
4=0,^-—1 показывает, что элемент х с одним и тем же индексом "- (1 -4) в: 


верхних левых углах блоков повторяется К раз. 
Для демонстрации шага индекса между секциями, шага индекса внутри секции 
и корректности построения свертки достаточно использовать более компактное 


представление: 

Ус № Хмь -.. | Хм, Хмрнк --. 

У к № | Ами Хм | 

о. й Я и г | И * с. о | 

У, х М № Хм 

У ХК Х, ЖЖ № ! 

| | 

1 Я 1 

| | 


[ 7, г! | ['; 7 з И 
Рисунок 3 — Индексы в секции и между секциями 


Из рис. 3 несложно определить количество элементов в каждой секции и число 
секций, необходимых для построения свертки. Из рис. 3 также видно, что с учетом 
того, что шаг индекса внутри каждой секции равен А и шаг индекса между секциями 
делится на 7, где г=А+1, свертка построена корректно. 

С учетом (12) и операторов И и 5 вычисление свертки М=К.-и, г=К-1, 
можно представить следующим образом: 


075, ох, 5, Ху, р 

Об. бу О бы О Жи. 
| 

ВЕ ВОЕН ЕН Е ИЕ ЗЕНОН ЕВА 

В 2х 0°5°Х, 

в. 


Рисунок 4 — Вычисление свертки длины М=К-и, г=А+1 в операторном виде 


> 7 Ге 7 Ге 7 ТС 
С учетом обозначений Ху =И'5,Ху, У =0'5,У,, Ку =И'5,Ку рис. 4 можно 
представить в блочном виде, в котором вычисление свертки длины №М=К.и, 
(К,’) =1, и-А =1, сводится к вычислению К сверток длины г. 


Лемма доказана. 
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Рисунок 5 — Свертка длины М№=А-и, Г=К +1 в блочном виде 


Продемонстрируем вычисление свертки при К - нечетном, г=А+1 и покажем 
на примере вычисления свертки длины № =12 =3.4. Вычисление можно представить 


в виде: 
У № № Ммм № № мм М № № 
( 5, У | М № ЖА № № ФМ № ЛМ 4% 
У № № № №№ М № №№ № м М 
а О о ое а И 
Ум м №№ Ммм № № №№ № № ХХ, 
сл Я у. Ум м м № М № ЖЖ м № № № 
Ув | Мл М м ММ м М № М № 
В Е № К И: ВЕ 
№ № № мм М № ММ № № № 
(7? 52, м № № Хм м М №№ М № Ж№ 
УХ м № №10 м м ми № М № 
У 1% № м № № № м № № № 
ое Ве а о М. о а 
05. В. 5, У В. 
Рисунок 6 — Вычисления свертки длины М№М=3.4 

Оо, И’ Ь |’ | И’ 

ПИ м. ОА М 

О о, АБ 09° Хр 

Вы: | ОВ | 5. В, 


Рисунок 7 — Вычисление свертки длины ЛМ =3.4 в операторном виде 


Индексы в операторах 5 $ и для входных и выходного сигналов совпадают в 


соответствующих секциях. 


С учетом того, что для вычисления свертки длины М№М=3 и № =4 необходимо 
4 и 5 операций умножения соответственно [6], для вычисления свертки длины 
№=12=3.4 необходимо 20 операций умножения, что соответствует оценке 
Винограда [8]. Это на 4 операции меньше, чем разбиение М =12=2.6=2.(2.3) на 


две секции по 6 элементов, которые разбиваются далее на 2 секции по 3 элемента 
(для вычисления свертки таким способ необходимо 24 операции однословного 


умножения [4]). 
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Лемма 5. Циклическую свертку А, длины М№М=КА.г, Г-К=А>1, (К,Г)=1 
сигналов Ху и У, соотношения (р.4), РА = ЙО А: РОТ 
РАЙ КЮ ре (р й. ‚ [= 0, К —1 переводят к блочно-циклическому виду. 


Доказательство. Воспользуемся оператором 5 (8) и найдем секции длиной г 
входных Ху, Ух и выходного Ку сигналов. С учетом леммы 2 элементы каждой 


секции сдвинем циклически вверх р, раз, 1 = 0, 


р 
2=1+р-й, р=а" =((-Ю7); р =(- р), д =(2),1=0Ж-Т; 


й 
Г 


№ к: к 
Хх, ХЕ Х + 
а Ро а Е, 
ОЕ Х:, 2% ЕО 
| Хм ы Х,, =й 211 
С учетом (13) вычисление свертки можно представить в виде: 
05°, пы Як ОХ, 
и А. | 0°5°Х, И (7 52Х, 
| сл _ 
... ... ы ... | | ... =. 
о о о 
ПА. ОА. Сы. О 


Рисунок 8 — Вычисление свертки с использованием операторов И и 5 


Вычисление свертки можно представить в виде следующих выражений: 
К 7— 


1 
—1 ы ь = г 
р (@ —Ю) й ‚ 2=1+р.К; Ка), = 1 88 енто, "Увы, 
= = 
РЕ (2 | ), › 
1=0, А-1, т=0,^-1, 4=0, -1. (14) 


Из (14) видно, что в каждой колонке (4 -7+т К) свертки индексы при хи у 


в каждой строчке (2, +7: К) „ отличаются на одно и то же число ({- а:г-т: К) ЧТО 
говорит о свойстве цикличности данного представления. 
я 1 _ТГРЕ © 7 _ТТР: © В р: <" 
С учетом обозначений Ху =И”5,.Ху, У =И"5,У,, В =И”5,Ву рис. 8 
можно представить в блочном виде, аналогичном рис. 5. 


Лемма доказана. 
Лемма 4 является частным случаем леммы 5 и (14) равнозначна (12) при 


г-К=а=1 (р=(4”"). =(-ю”). =1, 2=Р). 
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Априорные оценки сложности 


Теорема 1. Априорная оценка сложности вычисления циклической свертки 
длины М=К-к, (К,г) =1, "К =А>1, где К иг - взаимно простые числа, имеет вид: 


О`(№)< 0'(®)-0'(), 0*(№)< 9" (®-г+ 0-0“), (15) 
где О’(®), О"(+), О*(®, О*(+) - количество операций умножения, сложения 


и вычитания для вычисления циклических сверток длины К иг. 
Доказательство. Циклическая свертка сигналов Ху и У), где М=КА.г, 


(К,г) =1, "-К=а >1, может быть представлена в блочно-циклическом виде (рис. 5, 


8). Если О"(К), О"(ю), О*(® , О*(+) - количество операций умножения, сложения и 
вычитания, необходимых для вычисления циклических сверток длины К и г, 
то общую сложность вычисления свертки длины М№М=А.’ можно представить 
априорными оценками (15). 

Теорема доказана. 

Теорема 1 применима в следующих случаях: 

1. К — нечетное, г — четное, А <г. 

2. К — нечетное, г -— нечетное, К <г, К иг - взаимно простые. 

3. К — четное, г — нечетное, А <и. 

Случай, когда КА, г — четные, не рассматривается, так как по условию числа А и 
’ взаимно простые. Вычисление свертки вида М=2.К, К - нечетное [4], является 
частным случаем теоремы 1. 

Теорема 2. Априорная оценка сложности вычисления циклической свертки 


п-1 
длины № = | м,, (тт, =1, 12], т <т, <...<т, < т, у, имеет вид: 
1=0 


1-1 


ом) <П 9’, 0'м)< У | Пт) |.0*@т)-| Пт || 


К=0 К=+1 
* + < 
где О (т,), О”(т,;) - количество операций умножения, сложения и вычитания, 
необходимых для вычисления циклической свертки длины т, . 


Доказательство. Рассмотрим свертку длины М№М=т,.М№, где (ть,М)=1 


взаимно простые, тогда согласно теореме 1 априорные оценки сложности можно 
выразить следующими выражениями: 


О`(№<О’(т,)- 9°(М,), О" (М) < О* (т) М +0’(т)- 9" (м), (16) 
где О’ (т), О*(ть), О’ (№), О*(М№,) - количество операций умножения, 
сложения и вычитания, необходимых для вычисления циклических сверток длины 
ИА 
Пусть №=т ‘Л также является произведением взаимно простых чисел 
(т,,№,) =1, тогда согласно теореме 1 и формуле (16) априорные оценки сложности для 


свертки длины Л = т, - № = т т, . М можно выразить следующими выражениями: 
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0'(№)<О0"(т,)-0’(М)=0’(т)-(0"(т,)-9°(м,), 
0*(№) < О*(т,)-М, +0'(т,)- О" (М) 
0*(№) < О*(т,)-т,- М +0"(т,)-(0*(т,)- м, +0"(т)-0* (№) 
0*(№)<О*(т,)-т,- М +0"(т,)-О*(т,)-М, +0" (т,)-О"(т,)-9* (№), 
где О"(т,), О*(т,), О’(т,), О*(т)), О"(М,), О*(М№,) — количество операций 


умножения, сложения и вычитания, необходимых для вычисления циклических 
сверток длины т, т и М. Продолжая разбиение Л, на меньшие длины, 


приходим к соотношениям (16). 
Теорема доказана. 
Оценки в теоремах | и 2 аналогичны оценкам, полученным Р.К. Агарвалом и 


Дж.У. Кули при реализации метода с использованием китайской теоремы об 
остатках [9, с. 104-105]. 
К-1 
Следствие 1. Если сомножители являются степенью числа: № =[ [т ', То 
1=0 
деление на меньшие свертки должно производиться в порядке возрастания величины 
сомножителей т,’ <т,' <...<т,“? <т, 1, а не в порядке возрастания основания 


числа т <т <...<т_.<т, |. Так разбиение свертки длины 2”.3“.5”.7” про- 


изводится следующим образом: Вы (25<32<49<81). Сначала каждый 
сигнал разбивается на 25 секций, каждая из которых разбивается на 32 секции и т.д. 

Следствие 2. При группировке сомножителей порядок разбиения меняется. 

Имеется в виду, что число сомножителей может быть уменьшено за счет 
группировки сомножителей. Так, в предыдущем примере при группировке сомножителей 
25 и 49 меняется порядок разбиения 32<81<25.49. Это свойство является очень 
важным, так как для одной и той же длины свертки возможны различные разбиения с 
различной оценкой сложности, среди которых можно выделить лучшую оценку. Важным 
является то, что возможно комбинирование методов. Так, например, в предыдущем 
примере 32 <81<25.49 первый сомножитель является степенью двойки, что позволяет 
использовать методы, основанные на представлении большей свертки за счет нескольких 
сверток половинной длины. Для остальных шагов ...<81<25.49 будут использоваться 
другие методы. Если в начале последовательности 2<3<5<7<... стоят маленькие 
числа, то актуальным встает вопрос оптимизации вычислений не только длинных, но и 
коротких сверток. 


Таблица 6 — Количество операций умножения для вычисления коротких 
сверток [6, с. 67-70] 


т: ЕЙ ВЕН ВИ 
ом |245 0| 6 м |1 


Далее в таблицах приведены априорные оценки сложности по количеству 
операций однословного умножения для вычисления циклической свертки различной 
длины разными методами. 
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Таблица 7 — Количество операций умножения для вычисления сверток длиной, 
равной произведению взаимно простых чисел предложенным методом 


М = | [т, 6=2.3 | 10=2.5 | 18=2.9 | 24=3.8 За 7 63=7.9 


п-1 
О"(№) = ПО’(т) | 8=2-4 | 20=2.10 | 38=2.19 | 42=3.14 | 160=10-16 | 304=16.19 
1=0 


Таблица 8 — Количество операций умножения для вычисления сверток 
с использованием преобразования Уолша [3] 


п 3 4 ь 6 


М№М=2” 8 16 32 64 


О'(№=5.3"? | 15 | 45 135 | 405 


С учетом того, что при умножении двух М№ разрядных чисел получается 2№ 
разрядное число, можно использовать свертку длины 63 для умножения чисел 
длиной 32 слова. Согласно табл. 7 необходимо 304 целочисленных операций 
однословного умножения для реализации операции 32 разрядного умножения. Для 
сравнения: для вычисления дискретного преобразования Фурье длины 64 для одного 
числа необходимо (64/2)-105, 64 =192 операции, где присутствуют операции с 


плавающей запятой, а для реализации операции умножения, основанной на пре- 
образовании Фурье необходимо выполнять прямое и обратное преобразования. 

Метод удобен также тем, что пред- и поствычисления в виде циклических 
сдвигов для обоих входных сигналов одинаковые, что является отличием от методов 
Агарвала — Кули и Агарвала — Барраса. 


Выводы 


В данной работе предложен метод вычисления сверток длины, равной произве- 
дению взаимно простых чисел с использованием циклических сдвигов и без использова- 
ния китайской теоремы об остатках. Предложенный метод позволяет параллельно 
вычислять свертки меньшей размерности. Приведены априорные оценки сложности 
вычисления сверток длиной, равной произведению двух и более взаимно простых 
сомножителей. Приведены зависимости между элементами секций, числом секций, 
длиной секций, на которые разбиваются свертки. Приведены примеры разбиения цикли- 
ческих сверток разных длин. Предложенный метод может быть использован для вычи- 
сления циклической свертки целочисленных, вещественных или комплексных сигналов. 
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А.М. Терещенко, В.К. Задрака, А.М. Кун 
Використання цикл!чних зсув!в для прискореного обчислення циклёчних згорток довжини, 
яка дорвнюе множенню взаемно простих чисел 
п-1 
В робот! розглядаеться метод обчислення циклчних згорток довжиною № = ” ‚де (тт ,) =1, 
1=0 
1з2 } . Розглядаеться нове блоково-циклчне представлення згортки за рахунок використання циклЁчних 
зсув1в у кожному блоц1. При такому шдход! китайська теорема про залишки (КТЗ) не використовуеться. 
Метод е ефективним за рахунок використання нескладних перед- та шсляобчислень 13 застосуванням 
циклчних зсувв. Ейн дозволяе розпаралелювати обчислення. Отриман! апр!орн! ошнки складност! за 
юалькстю операщй множення. Наведен1 приклади обчислення згорток. 


А.М. Тегезйсйенко, Г.К. Гафгака, А.М. Кит 
Озто оЁ СусИе 5Ы $ Гог Еа$ Сотршайоп о? Сусйе Сопуошйоп о? Гепо® Едиа! 
0 Мш@рНсабоп оЁ Ваайуеу Риитез 
п-1 
Тве сотршавоп тео оЁ сусНс сопуоавоп оЁ мо земез еасН оР1епей оЁ № = П т; ‚ уве (т,,т,) =1, 
10 
Г ] 15 сопяаегед. Тре пе\ МосК-сусйс сопуоайоп ргезетаноп \й изо оЁ сусйс $$ ш еасв Ыоск 18 
ргорозеа. Те СЫтезе Ветатаег ТВеогет 15 поё зе ш 415 те@оа. ТВе тле#о4 15 еЙесйуе де ю изше зпире 
рге- ап розсотршанопз НКе сусИс $$. № аПо\у$ сотриайоп$ ш рагаПе]. ТБе патбег оЁ па @р|савоп$ т 
ргорозе4 плефоа 15 э1уеп. ТВе ехапар!ез оЁсопуо[аНоп сотрийоп$ аге ргезещед. 
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